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ABSTRACT 

The nature of the extended hard X-ray source XMMU J061804.3H-222732and its surroundings is investigated 
using XMM-Newton, Chandra, and Spitzer observations. This source is located in an interaction region of the 
IC 443 supernova remnant with a neighboring molecular cloud. The X-ray emission consists of a number of 
bright clumps embedded in an extended structured non-thermal X-ray nebula larger than 30" in size. Some 
clumps show evidence for line emission at ^ 1.9 keV and ^ 3.7 keV at the 99% confidence level. Large-scale 
diffuse radio emission of IC 443 passes over the source region, with an enhancement near the source. An IR 
source of about 14" x 7" size is prominent in the 24 fim, 70 /im, and 2.2 /im bands, adjacent to a putative Si 
K-shell X-ray line emission region. The observed IR/X-ray morphology and spectra are consistent with those 
expected for J/C-type shocks of different velocities driven by fragmented supernova ejecta colliding with the 
dense medium of a molecular cloud. The IR emission of the source detected by Spitzer can be attributed to both 
continuum emission from an HII region created by the ejecta fragment and line emission excited by shocks. 
This source region in IC 443 may be an example of a rather numerous population of hard X-ray /IR sources 
created by supernova explosions in the dense environment of star-forming regions. Alternative Galactic and 
extragalactic interpretations of the observed source are also discussed. 

Subject headings: ISM: individual (IC 443) — supernova remnants — X-rays: ISM 



L EVITRODUCTION 

The energy release and the ejection of nucleosynthesis 
products by supernovae (SNe) events are of great impor- 
tance for our understanding of the physics of the interstellar 
medium (ISM). The mixing of the ejected metals with the sur- 
rounding matter is of special interest when a SN occurs in a 
molecular cloud, which may cause further star-forming activ- 
ity. 

Optical and UV studies of the structure of SN remnants 
(SNRs) have revealed a complex metal composition of ejecta 
and the presence of isolated high-velocity ejecta fragments 
interacting with surrounding media. The most prominent 
manifestations of this phenomena are the fast moving knots 
observed in some young "oxygen-rich" SNRs, such as the 
Galactic SNRs Cas A (e.g., ChevaUer & Kirshner 1979; Fe- 
sen et al. 2006), Puppis A (Winkler & Kirshner 1985), 
G292.0H-1.8 (e.g. Winkler & Long 2006), and also N132D 
in the LMC and IE 0102.2-7219 in the SMC (e.g. Blair et al. 
2000). 

Ballistically moving ejecta fragments of SNRs can be con- 
sidered as a class of hard X-ray sources. The prototype was 
observed in the Vela SNR (Aschenbach, Egger, & Triimper 
1995; Miyata et al. 2001). A massive individual fragment 
moving supersonic ally through a molecular cloud can have a 
luminosity Lj > lO-'' ergs s"' in the 1-10 keV band, and is 
observable with XMM-Newton and Chandra at a few kpc dis- 
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tance (Bykov 2002, 2003). Its X-ray emission is expected to 
consist of two components. The first one is thermal X-ray 
emission from the hot shocked ambient gas behind the frag- 
ment bow shock, with a spectrum of an optically thin ther- 
mal plasma of an ISM-cloud abundance. The second emission 
component is nonthermal; the interaction of fast electrons ac- 
celerated at the fragment bow-shock with the fragment body 
produces a hard continuum as well as line emission (X-ray 
and IR), including the K-shell lines of Si, S, Ar, Ca, Fe, and 
other elements ejected by SN. Detection of the X-ray line 
emission would help distinguish an ejecta fragment from the 
other possible source of hard continuum emission associated 
with a SNR, namely, a pulsar wind nebula (PWN). 

A young SNR of an age of a few thousands years inter- 
acting with a molecular cloud can produce hundreds of X- 
ray sources associated with isolated ejecta fragments. They 
should be particularly numerous in starforming regions like 
those in the Galactic center region, where young core- 
collapsed supernovae in or near molecular clouds are expected 
to be present in abundance. The expected observational ap- 
pearance of isolated ejecta fragments in a molecular cloud 
differs from what is seen in the Vela SNR. Ejecta fragments 
interacting with a dense molecular cloud are slowed down and 
crushed, and they are generally more bright. We will argue 
here that the X-ray emission spectra of ejecta fragments in a 
molecular cloud may be dominated by hard non-thermal com- 
ponents, because a powerful but very soft thermal component 
could be heavily absorbed. On the other hand, the spectra 
of fast supernova ejecta fragments propagating in a tenuous 
plasma, as it is the case in the Vela SNR, would be long-lived, 
less luminous and dominated by thermal emission. 

The present paper focuses on IC 443. This is a SNR of a 
medium age, estimated by Chevalier (1999) to be ^ 30,000 
years, for which the number of X-ray sources from ejecta 
fragments should be much smaller than in a young SNR (pos- 
sibly, only a few). It is, however, the best and most reliable 
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laboratory to study this phenomenon since there are only very 
few examples of clearly established SNR-cloud interactions. 

IC 443 (G189.1+3.0) is an evolved SNR of about 45' size at 
a distance of 1.5 kpc (e.g. Fesen & Kirshner 1980). Radio ob- 
servations of IC 443 (e.g. Braun & Strom 1986; Green 1986; 
Leahy 2004) show two half-shells. This appearance is proba- 
bly due to interaction of the SNR with a molecular cloud that 
seems to separate the two half-shells. The molecular-cloud 
material has a torus-like structure (Cornett, Chin & Knapp 
1977; Burton et al. 1988; Troja, Bocchino, & Reale 2006), 
that can be interpreted as a sheet-like cloud first broken by 
the expanding pre-supernova wind and then by the SNR blast 
wave. Plenty of evidence for shock-excited molecules in this 
region has been found (e.g. DeNoyer 1979; Burton et al. 
1988; Dickman et al. 1992; Turner etal. 1992; van Dishoeck, 
Jansen, & Philhps 1993; Tauber et al. 1994; Richter, Graham, 
& Wright 1995; Cesarsky et al. 1999; Snell et al. 2005). The 
complex structure of the interaction region, with evidence for 
multiple dense clumps, is seen in 2MASS images (e.g. Rho 
et al. 2001). Three OH (1720 MHz) masers were found in 
IC 443 (Claussen et al. 1997; Hewitt et al. 2006, and refer- 
ences therein). 

Soft X-ray maps of IC 443 based on ROSAT data (Asaoka 
& Aschenbach 1994) and recent radio observations (Leahy 
2004) suggest that another SNR, G189.6H-3.3, is seen in the 
IC 443 field (see also the XMM-Newton study by Troja, 
Bocchino, & Reale 2006). This makes the multiwavelength 
observational picture even more complex to interpret. 

The field of IC 443 was observed in X-rays with HEAO 
1 (Petre et al. 1988), Ginga (Wang et al. 1992), ROSAT 
(Asaoka & Aschenbach 1994), ASCA (Keohane et al. 1997; 
Kawasaki et al. 2002), BeppoSAX (Preite-Martinez et al. 
2000; Bocchino & Bykov 2000), Chandra (Olbert et al. 2001; 
Bykov, Bocchino, & Pavlov 2005; Gaensler et al. 2006; 
Weisskopf et al. 2007), XMM-JVewtoii (Bocchino & Bykov 
2001, 2003; Troja, Bocchino, & Reale 2006), and RXTE 
(Sturner, Keohane, & Reimer 2004). 

The X-ray emission of IC 443 below 4 keV is dominated 
by a number of thermal components (e.g., Petre et al. 1988; 
Asaoka & Aschenbach 1994; Kawasaki et al. 2002; Troja, 
Bocchino, & Reale 2006). The thermal-emission morphology 
is center-filled, with soft emission filaments visible at ener- 
gies below 0.5 keV. A gradient of X-ray surface brightness 
at the SNR limb was found, as well as strong variations of 
absorbing column density A^h. which indicates the complex 
molecular-cloud environment of IC 443 in the southern part 
of the remnant (e.g. Asaoka & Aschenbach 1994). 

ASCA observations have established that the hard X-ray 
emission of IC 443 (above 4 keV) is dominated by localized 
sources in the southern part of the remnant (Keohane et al. 
1997). In XMM-Newton observations Bocchino & Bykov 
(2003; BB03 hereafter) found 12 sources with fluxes over 
10"''* ergs cm"^ s"' in the 2-10 keV band. Six of the de- 
tected sources are located in a relatively small, of 15' x 15' 
size, region projected onto the molecular cloud in the South- 
Eastem part of IC 443. BeppoSAX MECS observations 
(4-10 keV) showed two sources, ISAX J06 17. 1-1-2221 and 
ISAX J0618.0-I-2227, wifli evidence from flie BeppoSAX 
PDS for the presence of hard emission up to 100 keV for 
the former (Bocchino & Bykov 2000). Observations of this 
source by Chandra (Olbert et al. 2001; Gaensler et al. 2006; 
Weisskopf et al. 2007) and XMM-Newton (Bocchino & 
Bykov 2001) established its plerionic nature. Leahy (2004) 
argued that the pulsar that powers this plerion is associated 



with G189.6H-3.3 rather than IC 443. The nature of the sec- 
ond hard source - ISAX J0618.0H-2227 - remained unknown. 
This source, the brightest in the region (excluding the ple- 
rion), was resolved with XMM-Newton into two sources - 
the extended XMMU J061804.3H-222732 (of ~ 20" size) 
and file point-Uke XMMU J061806.4H-222832. We will call 
them Src 1 and Src 2 respectively (note that the sources 
were listed as Src 11 and Src 12 in BB03). The position of 
XMMU J061804.3H-222732 in the remnant is illustrated in 
Figures 1 and 2. A dedicated Chandra observation of Src 1 
has revealed a complex structure of a few bright clumps em- 
bedded in extended emission of > 20" size (Bykov, Bocchino, 
& Pavlov 2005; BBP05 hereafter). The brightest clumps are 
the extended Src la and the point-like Src lb. The appar- 
ent position of the source in a SNR - molecular cloud inter- 
action region naturally leads to SNR-related interpretations. 
The observed X-ray morphology of Src 1 and the spectra of its 
components are consistent with expectations for a SN ejecta 
fragment interacting with a dense ambient medium. Alterna- 
tively, Src 1 could be interpreted as a PWN associated with 
either IC 443 or G189.6h-3.3 (BBP05). However, one cannot 
exclude the extragalactic origin of the source, that is discussed 
in some detail in Section 4.7. 

IC 443 is a candidate counterpart of the EGRET 7-ray 
source 3EG J0617H-2238, wifli a flux of about 5 x 10"^ cm'^ 
s"' above 100 MeV (Esposito et al. 1996). The spec- 
trum of Src 1 extrapolated into the EGRET range is con- 
sistent with that of 3EG J0617H-2238 (BBP05). Also the 
position of Src 1 is consistent (albeit marginally) with that 
of 3EG J0617H-2238. Such a 7-ray luminosity can be ex- 
pected for both the fragment and PWN interpretations. The 
forthcoming GLAST mission (e.g. Johnson 2006) will 
be able to provide an accurate position and spectrum of 
3EG J0617H-2238, thus helping to solve the issue. Src 1 
Ues far away from the 99% error circle of the TeV-regime 
source recently reported by MAGIC (Albert et al. 2007) in 
the Western part of IC 443 field. The apparent position of TeV 
MAGIC source is close to the 1720 MHz OH maser detected 
by Claussen et al. (1997). 

We present here new results of a deep 80 ks observation of 
the region with XMM-Newton, imaging of the region with the 
Spitzer infrared observatory, and a new analysis of VLA radio 
observations. In § 2 a combined analysis of the new XMM- 
Newton observations and all the previous high-resolution X- 
ray data from Chandra and XMM-JVewton is presented, in- 
cluding images, spectra, and time variations in the X-ray do- 
main. In § 3 archival radio (VLA), IR (2MASS and Spitzer 
MIPS), and optical (POSS-11) data are used to constrain the 
nature of Src 1 . A discussion of the obtained results and fu- 
ture prospects are presented in § 4. 

2. X-RAY DATA ANALYSIS 

XMMU J061804.3H-222732 (Src 1) was observed wifli 
XMM-Newton in 2000 and 2006, and with Chandra in 2000 
and 2004 (Table 1). ClAO v.3.0 wifli CALDB v.3.2.2 was 
used for Chandra data processing, SAS v.20060628_1801- 
7.0.0 for XMM-Newton data processing, and the HEA- 
SOFT 6.1 suite, including XSPEC v.12.3, for spectral fitting. 

In the course of XMM-Newton data reduction, patterns 0-4 
of EPIC-PN detector events and patterns 0-12 of MOS events 
were used. To filter out the periods of enhanced particle back- 
ground, we applied a standard method of soft flare detection*. 

* http://heasarc.gsfc.nasa.gov/docs/xmm/abc/node7.html 
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Fig. 1 . — Wide-field views of SNR IC 443. Left: XMM-Newton 2-8 keV image with Spitzer MIPS 24 fim contours overlaid. Right: Spitzer MIPS image at 
24 fim with VLA 1.4 GHz contours overlaid. The images are produced from the XMM-Newton observations 01 14100101-01 14100601 and 0301960101, and 
Spitzer MIPS observations r4616960, r4617216, and r4617472. The white aiTow points to the studied region. 



The filtering did not significantly change the good time for the 
data taken in 2000, but reduced it in the dataset of 2006. 

In the Chandra observation of the year 2000, Src 1 is offset 
17. '6 from the ACIS aim point, on chip S4. To increase the 
quality of the ACIS chip S4 data the standard procedure was 
used to create an event2 file from the archival eventl file, in- 
cluding a run of the destreak task. As a result, the good time 
of the observation decreased from 11.5 to 9.6 ks. An appro- 
priate webscript^ was used to correct the Chandra event files 
for the systematic position offset. 

Figure 1 shows a wide-field view of IC 443, both in X-rays 
and in the medium IR 24 /^m band. Two X-ray images of 
the source region of interest here, obtained with Chandra and 
XMM-Newton , are presented in Figure 2. The coordinates of 
Src la and Src lb obtained from the Chandra data with the 
celldetect algorithm (a = 0648"04?32,(5 + 22°27'24f'l, and 
a = 06*'18'"04?58,(5 = +22°27'31"8, respectively) do not dif- 
fer from those reported by BBP05, as well as the position of 
Src 2 (a = 06^18™06?53,(5 = +22°28'28."2). A nearby source, 
which we refer to as Src 3, is located at a = 06''17™59?25, (5 = 
+22°27'38."9. In addition, the XMM-Newton data show ev- 
idence for an extended 'bridge' of diffuse emission between 
Src 1 and Src 3 (see Figures 1, 2, and 6). On the timescale 
of about 6 years, covered by the XMM-Newton observations 
used here (Table 1), it is not possible to detect a proper mo- 
tion of the sources at a 1.5 kpc distance if they move with a 
transverse velocity < 3000 km s"'. 

2.1. Source spectra 

The line of sight to the source region intersects the super- 
nova shell(s), the interior of the remnant, and the molecular 
cloud, each having different physical parameters. Therefore, 
the X-ray emission detected along the line of sight has mul- 
tiple thermal and nonthermal components, in addition to the 
Galactic background. A SNR shell can be characterized by 
a power-law emission spectrum, the hot low-density interiors 

' http://asc.harvard.edu/cal/ASPECT/fix_offset/fix_offset.cgi 



— by thermal plasma emission, and a dense molecular cloud 
may cause an appreciable absorption. 

We have performed a series of spectral studies of Srcs 1, 2, 
and 3, the results being presented in Table 2 (power-law mod- 
els). Table 3 (thermal models), and in the text of this subsec- 
tion. Mekal thermal plasma models were used (e.g. Kaastra 
1992) and the wabs model for absorption was applied (Mor- 
rison and McCammon 1983). Because of strong gradients in 
the soft X-ray surface brightness, one should carefully choose 
regions for spectral analysis. For the Chandra analysis, cir- 
cular source regions of 10" radius surrounded by 10" -wide 
background annuli were used. For XMM-Newton, whose 
PSF is broader, source regions of 20" radius and annuli with 
10" width were used. That results in an underestimation of 
a source flux caused by subtraction of the source photons re- 
garded as background photons due to the wide XMM-Newton 
PSF wings. The effect is estimated to be about 20% of the 
source flux. That may account for the differences in the Chan- 
dra and XMM-Newton normalizations of the fitted source flux 
reported in Table 2. However, the fitted power-law indices and 
absorption values are not affected. 

When comparing the Chandra and XMM-Newton spectra, 
such enlarged source extraction regions were also used for 
the Chandra data. For the spectral fitting, the count rate spec- 
tra were grouped with a minimum of 15 counts per bin. Ta- 
ble 2 summarizes the results of spectral fitting for Src 1, Src 2, 
and Src 3 with an absorbed power-law (using the wabs model 
for absorption; Morrison and McCammon 1983). Thermal 
plasma models for Src 1 and Src 3 yield very high and poorly 
constrained values of the temperature. However, the spectrum 
of Src 2 can be as well described with a mekal model as an ab- 
sorbed emission of thermal plasma (see Table 3). 

The high angular resolution of Chandra ACIS is useful for 
studying the spectra of the substructure of Src 1 in more de- 
tail. A 2" X A" elliptical source region around Src la, and 
a 2" radius source region around Src lb were selected. For 
both regions, the background counts were extracted from an 
annulus with inner radius of 10" and outer radius of 20" cen- 
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TABLE 1 

X-RAY OBSERVATIONS OF XMMU J06 1 804.3+222732. 



Obs ID 


Observatory 


Instrument 


Date of observation 


Exposure 


Good time 








(YYYY-MM-DD) 


(ks) 


(ks) 


760 


Chandra 


ACIS-S 


2000-04-10/11 


11.5 


9.6 






PN 




23.2 


19.5 


0114100301 


XMM 


MOSl 


2000-09-27 


25.6 


25.1 






MOS2 




25.6 


25.1 


4675 


Chandra 


ACIS-S 


2004-04-12/13 


58.4 


56.2 






PN 




79.9 


52.9 


0301960101 


XMM 


MOSl 


2006-03-30/31 


81.6 


67.6 






MOS2 




81.6 


68.2 



4Q.Q 

Src2 
30,0 ^ 

20.0 

10.0 



22:28:00.0 
50.0 

40.0 

4 

30.0 



20.0 

07.0 06.0 
27:10.0 



Chandra AC!S 

hardness ratio: 3 - 7 keV / 0.5 - 3 keV 




03.0 



02.0 



Src3 



01.0 6:18:00.0 17:59.0 



XMM 2-8keV 




22:00.0 

olu'mp D 

20:00.0 
■18:00.0 

. . 10 . iB ono.sno mo 



Fig. 2. — X-ray images of the region under consideration. The data are from the Chandra observation 4675 and from XMM-Newton observations 01 14100101- 
01 14100601 and 0301960101. Left: Chandra hardness ratio image (a 3.0-7.0 keV countrate map divided by a 0.5-3.0 keV countrate map) together with a contour 
outlining the extended near-IR source 2MASS J06180378+2227314 (with the L-shaped Si K-stiell X-ray emission line region inside). Broad band Chandra maps 
of the region have been presented by BPP05. Right. XMM-Newton image in the 2-8 keV band with superimposed contours of 2.122 /im H2 emission of Burton 
et al. (1987) [adopted from van Dishoeck, Jansen, & Phillips (1993)], indicating the presence of a shocked molecular cloud. The studied X-ray sources are 
indicated as well as the clump D of van Dishoeck, Jansen, & Phillips (1993). 



tered between the two sources. The spectrum of Src la can 
be described as an absorbed (Nh = (l.l±0.8)x 10^^ cm"^) 
power-law with photon index F = 1 .SilijJ ^ and a thermal plasma 
component with T = 0.2^g j keV. The reduced of this fit is 
0.71 at 14 d.o.f. The spectrum of Src lb can be described as 
an absorbed (Nh = 1.2;!;Q gX 10^^ cm"^) power-law with pho- 
ton index T = 2.0![j;^. The reduced of this fit is 0.73 at 9 
d.o.f. The quoted errors are at the 90% confidence level (for 
one interesting parameter). The confidence contours for the 
spectral parameters of the whole Src 1 are shown in Figure 3. 
These results are in a good agreement with those obtained by 
BBP05. 

The best-fit value of A^h for a power-law model of Src la 
(that dominates the 10" radius Src 1 in Table 2, line 4) is well 
below that for the South-Eastern region of IC 443 obtained 
with ROSAT and XMM-Newton (7 x 10^' cm'^). If one as- 
sumes that Src la, Src lb and the extended hard emission are 
at the same column density A^h'^ 7x10^' cm"^, one has to add 
a soft thermal component (of T > 0.1 keV) for Src la (See 
Fig. 3 in BBP05). Results of such fits, with the additional op- 
tically thin thermal plasma emission (mekal) component for 
Src la, are presented in Table 1 of BBP05. Such a component 
is localized in Src la and can be explained if the source is a 



supersonic ejecta fragment (given a range of shock velocities 
up to 300 km/s). 

Moreover, BBP05 showed that a satisfactory fit can be also 
obtained with a blackbody component used instead of mekal. 
Such a blackbody plus power-law fit, with a temperature of 
about 0.1 keV, radius of about 3 km, and power-law index 
of about 1.3 (at fixed A^h= 7x lO^' cm'^) can be understood 
— if Src la is a pulsar wind nebula — as a combination of 
thermal emission from the surface of a neutron star plus non- 
thermal emission from the neutron star magnetosphere and/or 
surrounding pulsar wind nebula. 

It should be noted that the high plasma temperatures ob- 
tained to fit a thermal model to the Chandra spectrum of the 
region of a 20" radius around Src 1 are due to a contribution of 
hard spectra of Src la, Src lb and their neighbourhood. If one 
considers the spectrum of this circle region with a removed 
inner circle of 1 1" radius surrounding Src la and Src lb (the 
resulting annulus is named Src 1* in Table 3), one obtains a 



1+0.5, 



X lO''^ cm-^ T = 5.3 



+11/ keV with x^ = 



fit with A^H= l-l-()4 
0.93 at 78 d.o.f. 

A line feature centered at « 1.8 keV is seen (Figure 4, top 
row) at the 99% confidence level in the XMM-Newton spec- 
trum extracted from the extended nebula of Src 1 — the L- 
shaped region shown in the left panel of Fig. 2. The combined 
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Fig . 3 . — Confidence contours (68%, 90% and 99% confidence levels) for the spectral fit parameters of Src 1 spectra obtained with Chandra. Upper and central: 
Temperature and photon index vs. hydrogen column density for a two-component thermal plasma + power law model. Lower: Photon index vs. hydrogen column 
density for an absorbed power law model. The outer contour on the central panel runs outside the shown parameter space. 



MOSl + PN + ACIS spectrum of the L-shaped region can 
be modeled as an absorbed A^h=(0.3±0. 2) x lO^^cm"^] power- 
law with r=l.l+0 2 and a possible Si hne at 1.9±0.1 keV of a 
fixed 0.1 keV width. 

To investigate the Chandra spectrum of the point-like Src 3, 
source counts were extracted from a small circle of 2" radius. 
The background counts were taken from an annulus with the 
inner radius of 3" and the outer radius of 10" . The spectrum 



is shown in Figure 4 (bottom row). It contains a feature at 
3.7 keV that is possibly due to an Ar emission line. There is, 
however, an alternative possible interpretation of the line as a 
redshifted Fe K line that assumes that Src 3 is extragalactic. 

The XMM-Newton spectrum of the faint bridge be- 
tween Src 1 and Src 3 can be modeled either as absorbed 
(A^H=0.5^S;:^xl022 cm-2) thermal plasma (T = 5.5^|;^ keV) 



emission with 



0.79 at 74 d.o.f. or as an absorbed 
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Fig. 4. — Upper Left: X-ray spectrum of the L-shaped region (shown in Fig. 2) derived from a combination of the MOSl, PN (obs. 2006), and ACIS (obs. 
2004) data. The model (histograms) corresponds to an absorbed [A'h=(0.3±0.2)x lO^^cm"^] power-law with r=l.l;';[J , and a possible Si hne at 1.9±0.1 keV of 
a fixed 0. 1 keV width. Upper Right: 68%, 90%, and 99% confidence contours for the parameters of the 1.8 keV Si K-shell line for the L-shaped region, shown in 
Fig. 2, obtained from combined data of the MOSl, PN (obs. 2006), and ACIS (obs. 2004) detectors. Lower Left: X-ray spectrum of Src 3 extracted from Chandra 



observations. The spectrum is modeled as an absorbed [Wh=0.6!J!-^ X 10 
width. Lower Right: 
confidence level. 

(A^H=0.6;|;ogX 10^^ cm"^) power-law of photon index F = 
2.0!^:^ with xl= 0.74 at 74 d.o.f. 

The 20" -aperture Chandra spectra of the studied sources 
Src 1 , Src 2 and Src 3 can be simuhaneously modeled by an 
absorbed (A^H=0.8;I;of x 10^^ cm~^) power-law of the photon 
index T = 1.8±0.2 with xl= 1-1 at 315 d.o.f. or by a ther- 
mal plasma (mekal) model with A^h=(0.6±0.1)x 10^- cm"^, 
T=12^] keV with xl= 1-1 at 315 d.o.f. 

Using all the data available, we found that Src 1 has not 
shown a significant time variation, its flux being consistent 
with that originally obtained by BB03. Some evidence (at the 
90% confidence level) is found for time variation in the un- 
absorbed flux of Srcs 2 and 3, as illustrated in Fig. 5. The 
unabsorbed flux of Src 2 in the 0.5-10 keV band decreased 
between 2000/09 and 2006/03 with 99% confidence. In the 2- 
10 keV band the flux increased between 2000/04 and 2000/09 
and later decreased with 90% confidence. The unabsorbed 
flux of Src 3 increased between 2000/09 and 2006/03 with 
90% confidence. More observations are needed to firmly con- 
clude on the issue. 

3. RADIO, IR AND OPTICAL DATA ANALYSIS 
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i (, , and a possible Ar line at 3.7±0. 1 keV of a 0.2_[j , i 
90%, and 99% confidence contours for the parameters of the 3.7 keV Ar K-shell line for Src 3. All the quoted eiTors are at the 90% 



] power-law with r=1.8^| 



3.1. VLA data analysis 

A radio image of IC 443 was obtained from archival VLA^ 
data obtained in 1997 at 1465 MHz from observations in 
the C and D arrays. The interferometric image has an an- 
gular resolution of 14" 1 x 13 "2. To recover information 
at all spatial frequencies, the synthesis data were combined 
with single dish data from the survey at 1408 MHz carried 
out with the 100 m MPIfR telescope (Reich, Reich & Fiirst 
1990). The final high fidelity image has an angular resolu- 
tion of 37."8 x 34"3, PA = 46.5°, and an average rms noise 
^0.1 mJy/beam. The total integrated flux density over the 
whole IC 443, 5 = (58 ± 2) Jy, is in a good agreement with the 
total integrated flux obtained from single dish observations 
(S ^ 60 Jy, Mufson et al. 1986), which assures the accuracy 
of flux density estimates over selected portions of the SNR. 

Figure 1 (right panel) shows the contours of the A 20 cm 
emission of SNR IC 443 superimposed on a 24 fim Spitzer 
map. The radio map for the source region of interest made 
from interferometric data only is shown in Figure 6 as con- 

The VLA of the NRAO is a facility of the NSF, operated under coopera- 
tive agreement by Associated Universities, Inc. 
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Fig. 6. — Left: Contours of radio surface brightness at A=2() cm superimposed on the XMM-Newton X-ray image in the 2-8 keV energy band. The angular 
resolution of the radio data is 14" 1 x 13"2. The radio contours are plotted at 1.2, 2.4, 3.6, 4.9, 6.1, 7.3, 8.0 and 8.6 mjy/beam. Right: Contours of 24/im 
emission detected by MIPS Spitzer superimposed on the XMM-Newton X-ray image in the 2-8 keV energy band. 



tours on a 2-8 keV XMM-Newton X-ray image. This fig- 
ure shows that Src 1 lies at the periphery of IC 443, far from 
the main SNR radio shell (that is situated in the North-East 
of the remnant), but near a localized radio excess. At the 
angular resolution and sensitivity of the present data, no ra- 
dio continuum source, either point-like or extended, could 
be associated with any of the X-ray sources. The local ra- 
dio flux density, obtained by integrating the radio emission 
over the region containing 97% of the XMM-Newton counts 
of XMMU J061804.3-H222732, is (60 ± 3) mJy. 

3.2. Spitzer MIPS imaging and photometry 

The field of IC 443 was the target of Spitzer MIPS scan 
observations r46 16960, r46 16960, and r46 16960 performed 
on 2005 November 9 (PI: G. Rieke). The Multiband Imag- 
ing Photometer for Spitzer (MIPS) aboard the Spitzer Space 
Telescope (Werner et al. 2004) is capable of imaging and 
photometry in broad medium-IR spectral bands centered at 
24 /xm, 70 /um, and 160 /xm, and low-resolution spectroscopy 
between 55 /im and 95 /im (Rieke et al. 2004). The 24 /im 
band covers the range of 21.3 /xm - 26.1 /itm, the lOfim band 
covers the range of 61.5 /xm - 80.5 /xm. Srcs 1, 2, and 3 were 



outside the field of view of the Spitzer IRAC near-IR camera. 

We used the standard MOPEX 030106 software (Makovoz 
et al. 2006) to construct mosaic images and extract point 
sources from the archival BCD-level data (pre-processed by 
the S 1 3 pipeline) according to the recipes in the Spitzer cook- 
books^ and the MIPS data handbook'". The first frames 
of each sequence were ignored. A total of 5940 individual 
frames were mosaiced for each of the MIPS bands. The net 
exposure of the mosaic maps is equal to 65-92 frames (2.62 
s each) for the 24/im band and 12-15 frames (3.15 s each) 
for the 70 /im band (different parts of the map were obtained 
with different effective exposures). Outlier detection was per- 
formed to exclude moving and solar system objects. A wide- 
field Spitzer MIPS image of the whole remnant is shown in 
Figure 1. 

Using the APEX software suite (Makovoz et al. 2006), 
we detected two point-like sources, at a = 06*'18™04?0, (5 = 
+22°27'23" and at a = Oe'^l 8^04.^0,(5 = +22°27'33" . The 
sources were detected only in the 24 /xm band with fluxes 
1.77±0.07mJy (S/N = 8.2) and 1.89±0.07mJy (S/N = 8.0), 

' http://ssc.spitzer.caltech.edu/documents/datademos/ 

http://ssc.spitzer.caltech.edU/mips/dh/mipsdatahandbook3.2.pdf 
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Fig. 7. — Src 1 environment as seen by Spitzer MIPS. The sources detected by Chandra are shown as crosses (3-8 keV) and a dashed rectangle (1.8-2.0 keV). 
The black contour denotes the extended near-IR source 2MASS J06180378+2227314. The dashed circles denote the beamsize of MIPS. 



respectively. The latter source coincides with the absorbed 
near-IR source 2MASS J06 180406+2227345 detected only in 
the K, band (7 > 17.9, //> 17.0, = 15.3 ±0.2). The dered- 
dened (2.02 - 2.30 /im) flux of the source ranges from 
0.64 to 1.21 mJy (considering the uncertainty of the extinc- 
tion value), that is (1.2-2.2)x lO"'-' ergs cm"^ s"'. 

According to the extinction maps of Schlegel, Finkbeiner & 
Davis (1998), the total Galactic absorption towards Src 1 cor- 
responds to Ay « 6.0 ±0.7. This sets an upper limit because 
the source is only ~1.5 kpc away. Moreover, the method of 
Schlegel et al. (1998) is known to overestimate the extinc- 
tion for dense regions by a factor of 1.3-1.5 (e.g., Arce & 
Goodman 1999). Thus, it is likely that Ay ^ 3^, consistent 
with the earlier estimates by van Dishoeck, Jansen, & Phillips 
(1993). 

There is another bright point-like source, 
2MASS J06 180359+2227227 = HST N8JT007783, pro- 
jected onto the area of Src 1. The source is clearly seen in 
the optical and near-IR bands (B = 16.3 ± 0.4, V = 14.5 ± 0.3, 
7 = 12.50 ± 0.03, i/ = 1 1 .88 ± 0.03, K,= 11.73 ± 0.02). The 
dereddened flux of the source ranges from 20.1 to 27.6 
mJy, that is (3.6-5.0) x 10"'^ ergscm"^s"'. Most likely, 
2MASS J06 180359+2227227 is a foreground star. 

An extended excess of IR emission is seen in the 24 /im 
mosaic map (left panel of Figure 7), coinciding with the ex- 
tended emission region detected with Chandra ACIS in the 
western part of Src 1 and with 2MASS J06 180378+22273 14, 
an extended (14" x 7" ) source of near-IR emission 
listed in the 2MASS XSC catalog with observed isophotal" 
Ks = 12.86±0.12. The excess is also seen in the 70/im band 
(right panel of Figure 7), its apparent size being comparable 
with the beam size. The aperture photometry estimates of 
the excess (with aperture corrections applied) are 1 1 A^l "^ mJy 
(90% err.) ±1.1 mJy (3(T error of the pipeline) for the 24 /im 
band and 840!^^|] mJy (90% eiT.) ± 170 mJy (3ct eiTor of 
the pipeline) for the 70/im band. The apertures of a 6" and 
8" radius were used for the 24 fim and 70 /im bands, respec- 
tively. These values correspond to (2.8!o:35 ± 0.28)x 10"'^ 
ergs cm-2 s"' for the 24/im band and (9.7!^;^f ± 1 .96) x 10"'^ 
ergs cm"^ s"' for the 70/im band. Notice, that the latter value 
is actually an upper limit. Depending on the extinction value, 

" http://www.ipac.caltech.edu/2mass/releases/allsky/doc/sec4_5e.html 



the dereddened flux of 2MASS J06 180378+22273 14 is 6.4- 
10.7 mJy, that is (1.2-1.9)x 10"'^ ergs cm"^ s"' in the 
band. 

The extended source is seen neither in the J and H 
bands of the 2MASS survey (Skrutskie et al. 2006) nor in 
the archival blue, red, and infrared images of the POSS- 
II survey (Reid et al. 1991). With the daophot pack- 
age (Stetson 1987) integrated into the NOAO IRAF soft- 
ware suite, the following upper limits were obtained for 
2MASS J06180378+2227314: 6.6x10"'^ ergs cm'^ s"' in 
the 2MASS J band, 2.9 xlO"'^ ergs cm'^ s"' in the 2MASS 
H band, 8.0x10"" ergs cm'^ s"' in the POSS-II blue band 
(3750-5500 A), 3.0 x lO''^ ergs cm'^ s"' in the POSS-II red 
band (5900-7100 A), and 3.0x10"" ergs cm"^ s"' in the 
POSS-II IR band (7350-8750 A). The limits assume the esti- 
mated extinction for the source (Ay = 6). 

The point-like X-ray source Src 3 is also seen as a weak 
IR source in all the bands of the 2MASS survey, with J = 
17.3 ± 0.3, // = 16.8 ± 0.6, K,= \6.2± 0.4 (see lower panels 
of Figure 8). 

4. DISCUSSION 

Below we discuss the multiwavelength data obtained for the 
~ 1 .'5 region that includes Src 1 and Src 3 and some implica- 
tions of the data analysis. We will refer to this extended source 
as J0618. 

The results of the multiwavelength data analysis presented 
above can be summarized as follows. The hard extended X- 
ray source Src 1 consists of 'clumps', both extended (Src la) 
and point-like (Src lb), which emit hard X-ray continuum. 
There are signatures of a thermal component of temperature 
^ 0.2 keV in the spectrum of Src la, and of X-ray line emit- 
ting clumps side by side with an extended IR source. The 
X-ray morphology of the extended source J0618 is complex. 
Src 3 is point-like according to the Chandra observations. The 
XMM-Newton image in Figure 2 shows evidence for an ex- 
tended bridge of diffuse emission connecting Src 1 with Src 3. 
Longer observations are needed to establish firmly the possi- 
ble presence of this faint bridge and of an X-ray halo around 
Src 3 (see e.g. Figure 6). 

The X-ray spectra of Srcs la and lb and Src 3 are character- 
ized by power-law components with photon indices F ~ 1.5. 
The spectrum of the Src 1 region shows a feature at about 1.8 
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Fig. 8. — Upper panels: Src 1 environment in the 2MASS near-IR bands. Tlie source regions and positions are marked in tlie same manner as in Figure 7. Tlie 
wliite X-mark denotes tlie position of 2MASS J06180359+2227227. Lower panels: Surroundings of Src 3 in the 2MASS near-IR bands. The 10" radius circle 
denotes the region used for spectral analysis of the X-ray data. The contour inside the circle denotes the position of Src 3 as seen by Chandra in the 0.3 - 10 keV 
band at the 2.5 X lO^^ cps/pix level. 



TABLE 2 

Parameters of absorbed power-law models of Srcs 1,2, and 3 



Source(radius)/Background 


Observatory, year 


Source counts 




r 


Norm^ 


* 

? 

xr./dof 


Src 1 (20" )/bkg2 
Src 1 (20" )/bkg2 
Src 1 (20" )/bkg2 
Src 1 (10" )/bkgl 


XMM 2000 A 
XMM 2006 A 
Chandra 2004 
Chandra 2004 


1064/505/507 
4574/1899/2150 
2973 
1314 


5+"-:' 

05+8:? 

0.7+"- 

2+8:' 


1 2+"-:' 
1 3+"'^ 

-0.04 

1 5+"'2 

1:0+8:^ 


2.8!'^ 10-' 
3.2+g-^10-5 
6.7+2:210-5 
1.4+"'^ 10-5 


1.13/124 
1.05/418 
0.84/151 
0.97/73 


Src 2 (20" )/bkg2 
Src 2 (20" )/bkg2 
Src 2 (20" )/bkg2 
Src 2 (10" )/bkgl 


XMM 2000 A 
XMM 2006 B 
Chandra 2004 
Chandra 2004 


955/418/470 
-/944/1003 

1829 

1270 


-Sr- 

1.0+":^ 
1.2+8:^ 
9+8:- 




'^■^-0.1 

2 7+0.5 

2.6+8:^ 

2 4+8:^ 


9.4+^5^10-' 
5.2!*'^ 10-5 
12.23210-5 
9.3+5'^ 10-5 


0.94/107 
0.88/105 
0.98/94 
1.05/66 


Src 3 (20" )/bkg2 
Src 3 (20" )/bkg2 
Src 3 (20" )/bkg2 
Src 3 (10" )/bkgl 


XMM 2000 B 
XMM 2006 A 
Chandra 2004 
Chandra 2004 


-/235/225 
2730/940/1080 
1251 
620 


5+5''' 
0.6+"* 
0.3+8:^ 
7+8:^ 


2 1+2:8 

1 6+"'* 
-Q.3 

1 5+<'* 
-fl.§ 

1 o+<).4 

^•^-0.3 


I.O+IV0-' 

1.0+1:210-5 
1.7+1:^10-5 


1.15/26 
1.10/234 

1.10/66 
1.36/34 



Note. — 'A' means PN, MOSl, MOS2 data combined, while 'B' means MOSl and M0S2 data combined; 'bkgl' denotes a background annulus 
with inner radius of 10" and outer radius of 20" , 'bkg2' denotes a background annulus with inner radius of 20" and outer radius of 30" . The source 
counts in the third column are given as PN/MOS1/MOS2 for XMM and ACIS for Chandra. All errors quoted in the table are at the 90% confidence level 
for one interesting parameter. 

"A/;/ is in units of 1022 cm-2. 

''Normalization parameter is the spectral flux in photons cm-2 ^-i j^gV-' at 1 keV. 

keV with a flux of a few times 10"^ ph cm"^ s"^ at the 90% ergs cm"^ s"' for the POSS-II infrared band (0.7-0.9 fim), 



confidence level, which might be attributed to a Si K-shell 
line. A feature at 3.7 keV was found in the X-ray spectrum of 
Src 3 at the 99% confidence level, which might be attributed 
to an Ar K-shell line, unless the line is a redshifted Fe K line 
of an extragalactic source. A firm detection of the lines from 
the localized clumps with XMM-Newton is hampered by the 
presence of strong hard non-thermal continuum emission. 

An extended source of IR emission was found at the North- 
western edge of Src 1 with a dereddened flux in the 24 /im 
band of SpitzerMIPS of about 3x10"'^ ergs cm"^ s"' and an 
upper limit for the 70/im band of about 10~" ergs cm"^ s"'. 
The near-IR flux of the source is about 2 x 10"'^ ergs cm"^ s"' 
in the 2MASS Kj band. The source is not seen in the J 
and H bands of the 2MASS survey, nor in the POSS-II op- 
tical bands. The upper limits are 3x10"'^ ergs cm"^ s"' for 
the H band, 7x10"'^ ergs cm"^ s"' for the J band, 3x10"'' 



3xl0-'2 ergs cm-2 s-i jj^g pQSS-II red band (0.6-0.7 
fim), and 8x10"" ergs cm'^ s"' for the POSS-II blue band 
(0.4-0.6 fim). These upper limits will be used below to con- 
strain the model we propose for J0618. 

Earlier studies of molecular emission from the extended 
region of apparent interaction of the IC 443 SNR with the 
neighboring molecular cloud have indicated the presence of 
emission from both fast and slow shocks (e.g. Burton et al. 
1988; Dickman et al. 1992; van Dishoeck et al. 1993, Snefl et 
al. 2005). Molecular clouds are known to have highly inho- 
mogeneous internal structure (e.g. Blitz 1993). A molecular 
cloud consists of numerous dense clumps with a rather small 
volume filling factor, embedded in an interclump matter of a 
modest density of 5-20 cm"^. The presence of a wide range 
of dense molecular emission clumps down to about 1" scale 
has been established in the cloud around Src 1 by Richter, 
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TABLE 3 

Parameters of absorbed plasma (mekal) models of Srcs 1, 2, and 3 



Source(radius)/Background 



Observatory, year 



Source counts 



T,keV 



Norm" 



Xj/dof 



n+D 

'-59.6 
75.1+"'** 

24 6+^^'^ 
-10.8 

79.8 



I.H+"\,10-* 



Src 1 (20" )/bkg2 
Src 1 (20" )/bkg2 
Src 1 (20" )/bkg2 
Src 1 (10" )/bkgl 



XMM 2000 A 
XMM 2006 A 
Chandra 2004 
Chandra 2004 



1064/505/507 
4574/1899/2150 
2973 
1314 



0.5!o, 
0.4^5- 
0.6^ 
0.3+ 



79.9+' 



1.1/124 
1.1/418 
0.85/151 
0.93/73 



1.2™-? 10^ 



Src 1* (20" )/bkg2 



Chandra 2004 



1781 



5.3^ 



-35*^ 
3:1^:^ 



0.93/78 



1.9!«|10^ 
8.6!|il0-5 
1.8!^|lO-^ 
1.5™1l0-4 



Src 2 (20" )/bkg2 
Src 2 (20" )/bkg2 
Src 2 (20" )/bkg2 
Src 2 (10" )/bkgl 



XMM 2000 A 
XMM 2006 B 
Chandra 2004 
Chandra 2004 



955/418/470 
-/944/1003 

1829 

1270 



i.o| 

0.7™' 



1.1/107 
0.98/105 
1.2/94 
1.4/66 



10^ 

10-5 

,-5 



Src 3 (20" )/bkg2 
Src 3 (20" )/bkg2 
Src 3 (20" )/bkg2 
Src 3 (10" )/bkgl 



XMM 2000 B 
XMM 2006 A 
Chandra 2004 
Chandra 2004 



-/235/225 
2730/940/1080 
1251 
620 



3+ 

o.4:*«- 

0.4: 

0.6- 



ii2_ 



7 2+ 
'■^-6.0 

36.7+22I 
o 1+71.8 

4 3+3-7 



1.9+ 
4.0+ 

3.9!^!,-^10 

°-0.9 



-0.9' 

4.8 10-5 



1.2/26 
1.2/234 

1.1/66 
1.3/34 



Note. — 'A' means PN, MOSl, MOS2 data combined, while 'B' means MOSl and M0S2 data combined; 'bkgl' denotes a background annulus 
with inner radius of 10" and outer radius of 20" , 'bkg2' denotes a background annulus with inner radius of 20" and outer radius of 30" . The source 
counts in the third column are given as PN/MOS1/MOS2 for XMM and ACTS for Chandra. AH errors quoted in the table are at the 90% confidence level 
for one interesting parameter. Src 1* is the annulus region of Src 1(20" ) with the removed 11" radius circle containing Srcs la and lb. 

"Nfj is in units of 1 0^^ cm"^ . 

Mekal model normalization parameter: K = 4^^2(14.^)2 / netindV, where D is the angular size distance to the source in cm, He and nu is the electron 
and hydrogen density in cnT^, integration is done over source volume. 



Graham, and Wright (1995). 

The apparent coincidence of the extended structured non- 
thermal X-ray continuum source J0618 with local excesses of 
IR emission can be understood in the framework of a model 
of interaction of a molecular cloud with a fast ballistically 
moving object. If the object was ejected by a SNR, it could 
be either a massive fast moving ejecta knot or a pulsar wind 
nebula. In both cases shock waves will be driven into the 
cold matter of the cloud. Such a scenario and the expected 
properties of the emitted IR and X-ray radiation will be dis- 
cussed below. A large power [above IQ-^'' erg s"'] released by 
a fragment of velocity about 300 km s"' being decelerated by 
a dense cloud is emitted as UV photons. The UV emission 
will create an HH region surrounding the fragment. The de- 
tected Spitzer emission and all of the IR/optical upper limits 
are consistent with the continuum emission of an HII region 
of a temperature about 10^ K and a number density above 100 
cm"^. However, to explain the detected near-IR flux of the 
source in the 2MASS K, band one should consider an addi- 
tional emission component most likely due to line emission of 
shocked molecular hydrogen. The radio, IR, and optical emis- 
sion of the Hll region surrounding J0618 will be addressed in 
Sections 4.2 and 4.3. 

4.1. X-ray emission from isolated fast ejecta fragments 

An important distinctive feature of J0618 is the presence of 
extended non-thermal hard X-ray emission. This can be ex- 
plained in the framework of a model of interaction of a mas- 
sive isolated ejecta fragment moving with a velocity above 
300 km s"' through a molecular cloud. In the model by Bykov 
(2002, 2003), X-ray emission from an ejecta fragment of 0.2 
- 0.3 pc size (~ 30" -40" at 1.5 kpc, the angular size of such 
a fragment in IC 443), interacting with a molecular cloud was 
estimated. The fragment mass was assumed to be > 10~^Mq , 
containing ~ 10~^Mq of Si. The "knot" traveled through an 
inter-clump medium of a number density ~ 100 cm"^ with a 
velocity of ^ 500 km s"' (the corresponding postshock tem- 
perature is ~ 0.3 keV). The model predicts a hard continuum 
emission with a photon index F < 1.5 and 1.8 keV Si line 



emission with a flux of ~ 10"*' ph cin"^ s"' . The K-shell lines 
are excited by both non-thermal and thermal electrons. 

In the case of J0618, a slightly more massive fragment 
highly structured due to interaction with the dense molecu- 
lar clump would be more realistic. A range of sub-fragment 
velocities around the mean value ^ 200-300 km s"' is consis- 
tent with the age of IC 443 of about 30,000 year advocated by 
Chevalier (1999). Such a fragment would provide line fluxes 
of a few times lO'^ ph cm"^ s'^ consistent with the esti- 
mated fluxes of the putative Si and Ar lines (Figure 4) from 
J0618. The effects of clumping of the metal-rich ejecta would 
result in an intermittent spatial structure of X-ray hue emis- 
sion from the source. It is worth to mention that the intrinsic 
absorbing column AA^h of a fragment could be substantial ( 
> lO'^^ cm"^), especially if the fragment contains metal-rich 
ejecta material. 

The relatively low velocity and the high absorbing column 
of an ejecta fragment in a dense molecular cloud makes its ob- 
servational appearance to be very different from that observed 
in Vela shrapnel A. The observed emission of Vela shrapnel 
A is strongly dominated by an optically thin thermal compo- 
nent (of T about 0.5 keV) of a shock heated plasma (Miyata 
et al. 2001). Contrary to the Vela case, a thermal emission of 
a thin plasma of temperature ~ 0. 1 keV is not the main com- 
ponent in the observed X-ray spectra of J0618 because of the 
high absorption. A hard non-thermal emission dominates the 
X-ray spectrum of J0618. 

4.2. Energetics of an X-ray emitting SN ejecta fragment 

The models of an isolated X-ray source involving the in- 
teraction of a supernova blast wave with dense ambient mat- 
ter suggest that the keV emission is due to bremsstrahlung 
of shock-accelerated electrons (Bykov et al. 2000; Bykov 
2002). Synchrotron X-ray emission would require electrons 
with TeV energies accelerated by shocks of speed well above 
1,000 km s"', which is hard to expect for a middle-aged SNR 
in a dense ambient medium, unless the moving object is a pul- 
sar wind nebula. 

The radiative efficiency of nonthermal bremsstrahlung is 
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known to be low at keV energies. This means that to 
produce a hard X-ray continuum at a rate eV, an electron 
of energy dissipates energy at a rate e'e via Coulomb 
losses in a medium of an average charge Z, and e'r/ee ~ 
6.3 X 10"^ Z{E^/m^c^) — see e.g. Akhiezer & Berestet- 
sky (1957). Therefore, to produce the X-ray emission at 
the observed level of about 6 x \Q^^d\^ ergs s~' by elec- 
trons accelerated to 50-100 keV, the dissipated power must 
be about e'e « dl ^Z~^{m(^c^ /E^) ergs s~^ The total dissi- 
pated power e^'in = lOe'eryl}. Here rj-x = ry/O.l < 1 is the ef- 
ficiency of electron acceleration by a shock, providing eyn ~ 
10^' d\^Z~^ rjZ\ ergs s ' .The electron energy E^ = 50 keV was 
assumed in the estimation. Since Z < 10 in an oxygen/silicon- 
rich gas an enhanced metallicity of an ejecta fragment (depen- 
dent on mixing of the ejecta with the ambient matter) could 
somewhat compensate the bremsstrahlung inefficiency pro- 
viding e'kin ~ lO^^c/j 5 ergs s"'. Note that the power required 
to produce radio-emitting relativistic electrons in clump D is 
also just above 10^^ g ergs s~^ 

The upstream ram pressure power dissipated at the for- 
ward shock of a ballistically moving fragment, e'sh « 8 x 

1035 

«a3 V2 '"20 ^1.5 ^rgs s \ is the source of gas heating and par- 
ticle acceleration. Here r2o is the shock radius of J0618 (mea- 
sured in 20" ). The mechanical power e'kin in the shock model 
must not exceed e'sh, resulting in the condition n.^T,v\ r^oZrj-i > 
12 to be fulfilled. Therefore, ejecta fragments of velocity 
V2 ~ 2 moving through a cloud of density n,^3 ~ 0.2 will have 
e'sh ~ ekin. if Z > 2. Thus, the observed patchy structure of the 
X-ray emission that is apparent in the Chandra ACTS images 
of J0618 could be attributed to structured metal clumps of the 
ejecta of IC 443. 

How can one directly detect or constrain the power dissi- 
pated by the shock in the model of IC 443 ejecta fragment bal- 
Ustically moving through a molecular cloud? The gas temper- 
ature behind the standard single-fluid strong MHD shock can 
be estimated as T w 1.4 ■ 10^ (K). For a shock with ener- 
getic particle acceleration efficiency > 10% that we consider 
here the postshock temperature should be reduced by a factor 
of about 1.2 to account for the effects of the energy flux car- 
ried away with energetic particles (see e.g. Bykov 2002). The 
postshock gas cooling distance estimated by Hartigan, Ray- 
mond, and Hartmann (1987) is about 1.8 xIO'^ V2^^ cm. 
The shock of a velocity 300 km s~' in the molecular cloud is 
radiative with the cooling layer angular size of 0.15" n~j. If 
the extended structure of Src la seen in Figure 2 of BBP05 is 
indeed due to the thermal emission of hot postshock gas then 
~ 0.03 would explain both the extension of about 5" and 
the luminosity of soft thermal component below 10^^ ergs s~' 
discussed in BBP05. Most of the shock power, however, is not 
in the soft X-ray emission of shock heated gas, but rather in 
UV-optical emission of the radiative shock dominated by UV 
lines of OVI at 1035 A, Ly^ and He II at 304 A(e.g. Hartigan, 
Raymond, and Hartmann 1987). 

At the same time the UV photons produced by the radia- 
tive bow shock of ejecta fragment with a luminosity Luv = 
10^^ erg i~'L36, will be absorbed and reprocessed, mostly to 
IR emission, in an expanding HII region surrounding the bow 
shock of J0618. Note here that the shape of the HII region 
could be different from that of the non-thermal X-ray nebu- 
lae, but it should overrun the bow shock of J0618. From the 
apparent position of J0618 in the molecular cloud and the es- 
timated fragment velocity one may conclude that it entered 
the molecular cloud about 1,000 years ago. Assuming that 



all of the UV photons produced by the radiative bow shock 
are absorbed in a spherical layer of an homogeneous ambient 
matter, one can obtain the estimation of the HII region radius 

of about6xlO"L3f n^l^^ cm (see e.g. Spitzer 1978). In fact, 
the geometry of J0618 is more complex than a spherical HII 
region because of the motion of the extended emitting frag- 
ment and because of the strong inhomogeneity of the ambient 
molecular cloud down to the arcsecond scale as was observed 
by Richter, Graham, and Wright (1995). In the interclump 
matter of density below 100 cm"-' the scale size of the HII re- 
gion would exceed the apparent X-ray size of J0618. On the 
other hand, the number density of the clump D located in a 
close vicinity of J0618 was estimated by van Dischoeck et 
al. (1993) to be above 10^ cm"^ and the size of the HII re- 
gion is consistent with the size of the 24 /im emission excess 
around clump D [of about (2-3)x 10~" ergs cm"^ s~'] appar- 
ent in Spitzer image presented in the right panel of Figure 6. 
Although very simplified, such a model allows us to estimate 
the radio, IR, and optical emission of the shock-produced HII 
region. 

4.3. Radio-IR-optical continuum and line emission of the 
HII region 

The emission of the HII region consists both of continuum 
emission of a thin thermal plasma of a kinetic temperature 
about 10,000 K and of a rich emission line spectrum (see 
e.g. Spitzer 1978). The specific appearance and the chemistry 
of the mostly neutral photodissociation region in molecular 
clouds were reviewed by HoUenbach and Tielens (1999). In 
the standard case the interstellar HII regions are powered by 
massive luminous stars. In the case of J0618 the source of 
the ionizing radiation is an extended fast moving bow shock 
and the ionized region is likely unsteady. That makes an 
accurate modelling of the system rather complicated, so we 
shall present here some approximate estimations of the ex- 
pected IR/optical fluxes from the HII region. The observed 
24 ixm Spitzer MIPS, 2MASS and 1.4 GHz radio emis- 
sion, as well as the upper limits provided by the other Spitzer, 
2MASS and POSS-Il observations of J0618 were modelled as 
continuum emission of a hot ionized plasma of a temperature 
ranging from 8,000 to 20,000 K. The modelled fluxes were 
corrected for interstellar extinction for a wide range of Ay 
values below 8, since Ay 8 roughly corresponds to A^h~ 
1.5x10^^ cm"^ - the maximal value allowed by the fits to 
the Chandra data presented in Fig. 3. Thus all the optical 
and IR measurements of J0618 were used to constrain both 
the kinetic temperature of the HII region and the parameter 
F = «g3 X rjo, where ne3 is the electron number density in the 
HII region measured in units of 10^ cm"^, and r2o is defined in 
Section 4.2. It was found that the measured 24 /xm flux of the 
source can be explained by a thermal continuum of a 10,000 
K temperature HII region and F ~ 0.03 with Ay > 6. Since 
the continuum emissivity is cx n^o, the apparent patchy image 
of 24 yum emission could be due to the clumpy structure of the 
ejecta fragment. 

The continuum emission from the HII region is consis- 
tent with the upper limits obtained from Spitzer, 2MASS and 
POSS-11 observations of J0618. The radio flux density from 
such an HII region is about 30 mJy. The flux density ob- 
tained by integrating the 1 .4 GHz radio emission over a larger 
(an arcminute scale size) region with VLA was about 60 mJy. 
Therefore the 1 .4 GHz radio flux is not in a conflict with the 
HII region model for the parameter Y ~ 0.03. We found also 
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that the flux of the hydrogen recombination Hne from the 
HII region is below the upper Hmit of 3 x 10"'^ ergs cm"^ s"' 
obtained from the POSS-II red band (0.6-0.7 /xm) observation 
discussed above. 
The estimated parameter Y ~ 0.03 corresponds roughly to 
~ 160 cm~^. The bow shock radius estimated from X-ray 
image is r2o ^ 1 implying that the derived electron density in 
the bright part of the HII region in J0618 is dominated by the 
compressed postshock flow. Then the preshock ambient den- 
sity in the radiative shock can be estimated from «a < ne/4, 
providing the ambient gas density < 40 cm"-', generally con- 
sistent with that expected in the interclump matter of a molec- 
ular cloud. A substantial part of the powerful UV emission 
from the radiative shock of J0618 will irradiate the nearby 
dense molecular clump D providing another HII region of a 
surface brightness and high luminosity well in excess of that 
from J0618 that is clearly seen in Fig. 6. In the frame of the 
ejecta fragment model of J0618 we estimated that at least a 
substantial amount of 24 /xm flux from clump D [that is about 
(2-3) X 10"" ergs cm'^ s'M can be attributed to the HII region 
excited by the UV emission from the nearby source J0618. 

That model provides, however, only about one third of the 
observed 2MASS band flux of J0618, and that can not 
be simply relaxed with an appropriate HII region parameters 
choice. Thus, an extra contribution in the 2MASS band 
at the level of ~ 1.5x10"'^ ergs cm~^ s~^ is required. That 
contribution may come, most likely, from the emission lines 
of the shocked molecular hydrogen in J0618 and would re- 
quire the presence of C-type molecular shocks of velocity 
^ 30 km s~' in the close vicinity of J0618. Moreover, the 
atomic fine-structure lines could contribute to the 24 fim emis- 
sion detected by Spitzer MIPS if supernova ejecta drive J-type 
shocks of about 100 km s~' (and faster) into dense molecular 
clumps. We discussed above only the emission produced by 
the forward bow shock of the fragment, however slower re- 
verse shock will also be present and since the ejecta fragment 
body is likely very inhomogeneous the real structure is likely 
even more complex. Relevant line emission models will be 
briefly reviewed in the next subsection. 

4.4. IR-line emission in radiative shock models 

Atomic fine structure lines of [OI] (63 fim) and [Fell] (26 
/im) could dominate the emission in the Spitzer MIPS 70 /im 
and 24/im bands. The lines are known to trace fast radia- 
tive shocks in molecular cloud material. A comprehensive 
study of radiative shocks in interstellar clouds has been done 
by HoUenbach and McKee (1989; HM89 hereafter). They 
have shown that the intensity ratio of [OI](63 ^m)/[FeII](26 
^m) is about 10 for a shock wave in a cloud of a density 
«ci =10^ cm"^. They estimated the [OI] (63 /um) IR line inten- 
sity as /63 » 1.5 X 10"^ ergs s~' sr~^ for a shock of a velocity 
Vsh =150 km s"' and demonstrated that it scales roughly lin- 
early with HciVsh- The combination of optical line emission 
of OI at 6300 A, CII (2326 A), and Oil (3726 A) provides 
the dominating gas coolant in transparent systems of na < 
l(f cm~^. The line intensities were obtained by HM89 under 
the assumption of standard solar composition with account 
for the interstellar gas depletion to dust grains. The IR-optical 
line intensities from shocks driven by (and into) metal-rich 
ejecta could be produced with lower pre-shock number densi- 
ties Uci than those modeled by HM89, and should be re-scaled 
respectively. This can reduce the pre-shock density required 
to match the observed fluxes. 

Applying the HM89 model intensities to an IR source of 



an angular area A km (measured in units of 100 arcsec^) as 
discussed in section 3.2, one obtains a [Fell] (26 fim) line 
flux F26 ~ 2.4 X 10"^^Aioon3V2 ergs cm~^ s"^ Here ns is 
the pre-shock number density in units of 10^ cm~^ and V2 
is the shock velocity in 10^ km s~'. To reach the 24 /zm 
flux of about 3 x 10"'^ ergs cm"-^ s~^ as estimated from the 
Spitzer MIPS data, one would need n3V2Aioo "^1- As the 
estimated emission area is Aioo ~ 2, one needs nj,v2 ~ 0.5 
to account for the observed IR emission. The [Fell] (26 ^m) 
line could provide a sizeable part of the 24 /im Spitzer MIPS 
flux. The associated model flux of the [01] (63 /vm) line, 
^63 ~ 2.4 x 10~'^Aio() • n3V2 ergs cm"^ s"' is consistent with 
the upper hmit derived above. It is apparent from Figure 2 
that Src 1 is located near the shocked molecular clump D 
presented in a map by van Dishoeck, Jansen, & PhiUips 
(1993). Typical densities of molecular clumps, studied by 
van Dishoeck et al. (1993) are > 10"^ cm"^. The condition 
niV2 ~ 0.5 can be matched by either a slow shock of veloc- 
ity V2 ~ 0.3 in the outskirts of clump D, or by a fast shock of 
V2 ~ 1 in the interclump medium. 

Near-IR metastable lines of Fell at 1 .3 /xm and 1 .7 /xm have 
a rather flat shock velocity dependence above 100 km s~' . The 
intensity of the brighter Fell line at 1 .3 /xm scales roughly lin- 
early with nci, and it is about 0.2/63 . This estimate is close 
to, but still consistent with the upper hmit of the flux in the 
2MASS J band discussed in Section 3.2. A near-IR NI (1.04 
/im) line flux was predicted by HM89 to be at the level com- 
parable with that of [OI] (63 /im), and a SI (1.1 /im) line flux 
similar to that of Fell at 1 .3 /im. 

From the apparent lack of extended emission in the POSS-II 
red band (0.6-0.7 /xm) the flux upper limit for OI (6300 A) can 
be estimated as 3 x 10"'^ ergs cm"^ s"' (assuming the reason- 
able extinction to Src 1 as Ay « 6). The flux limit is below the 
value predicted by HM89 for the 01(6300 A)/[Fen](26 /im) 
ratio for a transparent system with shock velocity V2 > 1. 
However, the line ratio is increasing with shock velocity for 
V2 < 1, and, for example, a fast J-type (see Draine 1980 for a 
definition of J- and C-type shocks) radiative shock of V2 ~ 0.6 
in a medium of n3 ~ 1, could explain both the observed [Fell] 
26 /im line flux and the derived upper limit for the OI (6300 
A) line flux. Note that optical line fluxes from a radiative 
shock propagating into oxygen-rich ejecta material could be 
reduced by heat-conduction effects (e.g. Borkowski & ShuU 
1990). 

An important diagnostic IR-line ratio sensitive to the shock 
velocity is [NeII](12.8 /im)/[FeII](26 /im). The ESO VLT 
Spectrometer and Imager for the Mid-Infrared (VISIR) is an 
optimal instrument for observations in the two mid-infrared 
atmospheric windows: the 8-13 /xm N band with a 19'.'2 x 
19."2 FOV and the 16.5-24.5 fim Q band with a 32"3 x 32."3 
FOV (Lagage et al. 2006). The spectrometer could detect 
[Nell] (12. 8 /im) predicted by the radiative shock model of 
HM89 and the continuum emission from Src 1 HII region, 
thus discriminating between the two models. 

The model of a single fast radiative J-type shock de- 
veloped by HM89 underpredicts the near-IR flux (above 
10"'^ ergs cm"^ s"') observed in the band from the ex- 
tended source 2MASS J06180378H-2227314, closely associ- 
ated with the Spitzer MIPS 24 /xm excess (Figure 7). Thus, a 
combination of fast and slow shocks is needed to explain the 
IR emission from J0618. The same conclusion has been made 
in most of the studies of molecular emission of the south- 
western region of IC 443. Burton et al. (1990) argued that 



X-ray-IR sources in IC 443 



13 



both fast dissociative J-type and slow C-type shocks must be 
present in the extended southern cloud of IC 443 to explain 
the observed line fluxes [see also Snell et al. (2005) for a 
recent discussion]. Analyzing different shock models, includ- 
ing a time-dependent one, Snell et al. (2005) concluded that 
no single-shock model can explain the existing observations, 
and a range of shock velocities is required. 

4.5. Shocks in a molecular cloud driven by ejecta fragments 

A wide range of shock velocities in a molecular cloud is re- 
quired to explain the IR and X-ray observations of Src 1 . Slow 
molecular C-type shocks are required to explain the observed 
2MASS band emission. 

Molecular emission occurs when a fast ejecta fragment col- 
lides with molecular clumps of different densities, pro- 

-1/2 

ducing multiple shocks of velocities Vj oc p^^ . This can ex- 
plain the shock velocity range of 10 km s~^ < Vs < 100 km s~\ 
required to explain both the fluxes detected by Spitzer MIPS 
and 2MASS K, and the upper limits on the optical lines. 

In the considered case the presence of a fast shock of veloc- 
ity > 300 km s~^ seems to be required to explain the X-ray 
data discussed above. Moreover, it is apparent from Figure 1 
that while the 1.4 GHz radio emission is dim in the south- 
eastern part of the shocked molecular cloud in comparison 
with that at the north-east, there is a localized excess in the 
vicinity of molecular clump D adjacent to J0618. If the radio 
emission is the synchrotron emission produced by relativistic 
electrons in the GeV regime in a likely enhanced magnetic 
field of the dense clump D, then the model requires the pres- 
ence of a fast shock of velocity well above 100 km s~' to 
accelerate the relativistic particles in the vicinity of Src 1. 

The shock ram pressure estimations for multiple shock 
models were discussed by many authors (e.g. Burton et al. 
1990; Richter, Graham, and Wright 1995; Cesarsky et al. 
1999; Snefl et al. 2005). A fast shock of velocity ~50- 
100 km s"' incident on a clump of a moderate number density 
lO^-lO"* cm"^ would have a high ram pressure nv^ ~ 10^- 
10^ cm"-' [ km s"']^. It seems to be uneasy to attribute the 
estimated range of ram pressures to a uniform SN blast wave 
colliding with a cloud. A SN blast-wave ram pressure, as es- 
timated from the properties of shocked X-ray emitting gas 
inside the remnant, is smaller. On the other hand, strong 
ram-pressure inhomogeneities can be produced if a molecu- 
lar cloud is hit by SN ejecta driving the blast wave. The ejecta 
are likely to be rather fragmented. Each can be considered 
as an ensemble of ejecta fragments with a range of veloci- 
ties and densities produced by early ejecta instabilities. The 
fragments could provide a wide range of ram pressures in the 
molecular cloud (up to 10^ cm~^ [ km s~']^) driving shocks of 
different velocities, which are required to explain the observa- 
tional data. Isolated fast moving dense ejecta fragments could 
penetrate deeper into the dense parts of the cloud driving the 
high-pressure shocks on relatively small scales. 

4.6. Time variability and ejecta fragments statistics 

The Ufetime of an isolated ejecta fragment in a dense 
medium is an important factor with regard to time variabil- 
ity of the X-ray emission. A fast moving knot is deceler- 
ating due to the interaction with the ambient gas. The drag 
deceleration time of a fragment of velocity v, mass A4 and ra- 
dius TZ can be estimated as « 100A^-2(na3 V2 7?.^i )"' years. 
Here Tl-i = 7^/(0.1pc) and M-i = M/IO'^Mq . The num- 
ber density n^s of the ambient matter is in units of 10^ cm~^. 



and the fragment velocity V2 in units of 10^ km s Hy- 
drodynamical crushing of a fast knot occurs on a timescale 
Tc ~ x'/^7^v-' (e.g. ChevaUer 1975; Suflierland & Dopita 
1995; Wang & Chevalier 2001, and references therein). The 
density contrast x = Pk/ Pa (Pk is the knot density) is of the 
order of 1 for a large enough fragment, and ra < Tc for large 
enough clumps. To reach its apparent position in the molecu- 
lar cloud of IC 443, an isolated ejecta fragment must be mas- 
sive enough, A4-2 ^ 1, to overcome strong drag deceleration 
in the dense matter, as discussed below. 

Particle acceleration occurs on a few-years timescale if the 
fast particle diffusion coefficient does not exceed 10^^ cm^ s~^ 
(see Bykov 2002 for a discussion). Thus, given the fragment 
lifetime estimated above, some time variability of the X-ray 
emission, both in hard continuum and lines, can be expected 
on a few-years timescale (and longer) for a fragment of ve- 
locity V2 > 3 in a dense ambient matter with «a3 5^ 3, typical 
for molecular clumps. Variable X-ray emission has been re- 
cently observed with Chandra in a few molecular clouds in 
the Galactic Center region by Muno et al. (2007). Using the 
above estimates, one can argue that variable X-ray emission 
could be produced due to an interaction of metal-rich super- 
nova ejecta with dense clumps of molecular clouds. 

According to the log A'^ - log S distribution of X-ray emitting 
knots simulated by the method of Bykov (2003), the proba- 
bility of another similarly bright ejecta fragment getting into 
the field is rather low because of the short lifetime of the frag- 
ment in a dense molecular clump. However, some smaller and 
less massive fragments propagating in the interclump medium 
could be seen in the cloud as weak point-Uke X-ray sources of 
luminosities L^^ < 10^' ergs s~' . In order to firmly identify the 
sources as such fragments, a deeper observation is required. 

The hydrodynamic simulations available (Klein, McKee, & 
Colella 1994; Wang & Chevalier 2001) predict a complex ir- 
regular structure of the fragment body due to hydrodynamic 
instabilities. In this case the X-ray image would show an ir- 
regular patchy structure, instead of a smooth regular head-tail 
structure. The patches are due to emission of dense pieces 
of the fragmented knot illuminated by the shock-accelerated 
energetic particles. The observed morphology of X-ray emis- 
sion with bright clumps of a few arcsec size in Src 1, such as 
Src la and Src lb, as well as Src 3, can be explained in that 
model. It is worth to note that ballistically moving ejecta frag- 
ments could have rather large non-radial velocities (with re- 
spect to the apparent center of SN explosion). The non-radial 
velocity components will be substantial if the fragments origi- 
nate from ejecta instabilities at relatively late evolution stages. 

Thus, one can conclude that the morphology, spectra, and 
X-ray luminosity of the J0618 complex are generally consis- 
tent with those expected for a balUstically moving SN ejecta 
fragment interacting with a molecular cloud, though it is 
worth to discuss also some alternatives. 

4.7. Alternative interpretations of J0618 

A possible interpretation of the extended hard X-ray source 
J0618 associating it with an interaction of the ejecta of SNR 
IC 443 with the nearby molecular cloud has been discussed 
above in some detail. Another SNR-related possible interpre- 
tation is a low luminosity pulsar wind nebula originating ei- 
ther from IC 443 or from anoflier SNR G189.6H-3.3, advanced 
byBBP05. 

Extragalactic sources could also contribute to the observed 
appearance of J0618. Spectra the of point- like sources Src 2 
and Src 3 could be interpreted as that of AGNs. Within such 
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an interpretation the feature at about 3.7 keV in the spectrum 
of Src 3 could correspond to a redshifted (with 0.6 < z < 1) 
Fe emission Une. Moreover, one could speculate on the asso- 
ciation of the extended Src 1 with a high-z cluster of galaxies. 
The angular size of Src 1 (< 1' at photon energies > 2 keV) 
would suggest in that case a redshift z > 0.5 (see e.g. Rosati 
2004 for a review of high-z clusters). The thermal fit of the 
diffuse Src 1 * presented in Table 3 with the temperature about 
5 keV is consistent with the cluster interpretation. 

Optical and IR identifications of galaxies (of a few arcsec- 
ond scale size at z ^ 1) in the hypothetical cluster are the most 
obvious test to check the interpretation and measure the red- 
shift of the putative cluster. With the optical data available 
and given the substantial absorption (3 < Ay < 6), it is possi- 
ble to identify the objects not fainter than my = 17-19, while 
the brightest galaxies in most of the known X-ray clusters at 
z < 1 are of my > 19-20. Thus a dedicated optical-lR study 
is needed to confirm or reject the high-z galaxy cluster inter- 
pretation. The observed extension of the 2MASS Kj source 
correlated with observed diffuse X-ray emission is not easy to 
understand in the cluster of galaxies interpretation; it rather 
supports a Galactic SNR-related origin of J0618. 

5. CONCLUSIONS 

The multi-wavelength observations presented here indicate 
a possible physical connection of the X-ray source J061 8 with 
the neighboring IR Spitzer and 2MASS sources. That connec- 
tion, if real, can be understood in a scenario where J0618 orig- 
inates in an interaction of the IC 443 SNR with the adjacent 
molecular cloud. The correlation would require the presence 
of both fast and slow shocks in the clumpy molecular-cloud 
material. The X-ray line features apparent in the spectra of 
the clumps favor a scenario in which the shocks are produced 
by a fast ballistically moving SN ejecta fragment penetrat- 
ing into a structured molecular cloud. The model provides a 
physical picture coherent with the current observational data, 
although alternative scenario cannot be rejected yet. 

Alternatively, Src 1 can be interpreted as a massive X-ray 
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cluster of galaxies at a redshift z > 0.5. High-resolution 
arcsecond-scale observations and fine spectroscopy are re- 
quired to distinguish between these very different scenarios. 

If the SNR-ejecta interpretation is confirmed by further ob- 
servations, the source J0618 in the IC 443 could be a pro- 
totype of a rather numerous population of hard X-ray - IR 
sources created by SN explosions in the dense environments 
of star-forming regions. Such sources would be particularly 
abundant in the Galactic Centre region. 
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